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When noise radiation of roughly one or a few kilocycles per second is emitted in the 
higher atmosphere, part of it (an extraordinary wave) is propasjated downward into the 
space between earth and ionosphere. Reflection at the earth and ionosphere leads then to 
a standing wave in this space, whose intensity for a given incident powc^r flux varies very 
much with frequency. Maximum field strength of the standing wave is derived for the 
resonance frequencies of the space. The incident wave fronts are assumed to be horizontal. 
Only clearly defined wave fronts being of a sufficiently wide extension and showing no 
noticeable irregularities will lead to a definite resonance. The space between earth and 
ionosphere is comparable to an air gap between two parallel plane reflectors. The lower 
ionosphere is however only a partial reflector, allowing radiation to enter the resonance 
space and causing at the same time some loss of energy out of the resonance spacer (leakage). 
Stratification of the lower ionosphere has some* influ(MU'e on the resonance phenomena. 
With a D layer, an additional resonance at a lower fre(iuency is obtained. 

It seems too early to decide whether in any observed noises the resonance spectrum of 
the space between earth and ionosphere becomes apparent, but it is expected that noise 
spectra observed on the ground are modified by the resonances. In case of monochromatic 
emissions, the received intensity depends on the position in the resonance spectrum. The 
resonance spectrum should be received in case of emission of a white-noise spectrum, provided 
the wave fronts are appro})rijite. 



1. Introduction 
1.1. Cavity Between T\A<'0 Plates 

In a space between two infinite parallel plates 
resonance is o})served when half a wavelength or 
mnltiples of a half-wavelent>;th fit into the space. 
The electromagnetic field in this case of resonance is 
that of a standing wave whose wave planes are 
parallel to the plates. If one of the plates is made 
partially transparent, a model of the earth and 
ionosphere is obtained. The semitransparent plate 
like the ionosphere lets waves enter from ontside — or 
from the higher atmosphere — and still permits 
resonance oscillations. The phenomena enconntered 
in the case of two semitransparent plates (wire 
gratings) or one conductor together wath a semi- 
transparent plate (grating or dielectric plate) have 
been described in the literatm'e [Casey and Lewis, 
1952; Wait, 1954; von Trentini, 1956]. 

The field intensity in the resonance space is 
strongly dependent on the frequency of the incident 
waves, showing peaks at the resonance frequencies. 
The partial transparency of one of the reflectors 
entails imperfect reflection and contributes thus to 
a limitation of the resonance peaks. 

1.2. Fundamentals 

The ground is a good reflector for frequencies of 
the order of 1 kc/s. With a conductivity of 10~^ 
mho/m, for example, the amount of the complex 
refractive index of the ground for 1 kc/s is |7i| = 134, 
a value sufficient for good (almost perfect) reflection 
in vertical incidence, which is considered exclusively. 



The ionosphere is only a partial reflector. TIk^ 
extraordinai'v wave at very (or extr(Mnely) low fjc- 
quencies finds in the ionosphere (except its lowest 
part) conditions of free propagation, as is known 
from the audiofrequeiu'y whistlers. The refractive 
index of the extraordinary wave under these condi- 
tions is in general high (::^1). When the increase 
to these high values is fairly rapid in terms of wave- 
lengths, we may expect that part of the extraordinary 
wave passes through the height int(M'val of this rapid 
increase and part is reflected there. Thus wh^ have 
for extraordinary waves the combhiation of an almost 
perfect reflector (the ground) with a semitransparent 
reflector that allows waves to enter the I'csonance 
space. 

The descending wave fronts and consequently also 
the wave fronts in the resonance space will be as- 
sumed as horizontal. An eventual slight inclination 
of the wave fronts in the ionospliere w^ould, because 
of the refraction in leaving the ionosphere, cause a 
considcT-able inclination of the wave fronts in the 
resonance space with the result of an essential 
modification of the resonance spectrum. Random 
distortion of the wave fronts would correspond to a 
mixture of propagation directions and consequently 
lead to a smearing-out of the resonance spectrum. 
The lateral extension of the wave fronts must be 
wide enough to have the simple plane-wave theory 
applicable. 

There is some relationship between the resonance 
effect under consideration and the ^^tweeks", se- 
quences of pulses caused by multiple reflection of 
sferics in propagation to a distant location [Potter, 
1951]. The lightning flash from which a sferic 
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originates is nearly a point source. The radiation 
thus is omnidirectional and the tweek spectrum 
may not be interpreted as a resonance spectrum 
corresponding to horizontal (or otherwise simply 
defined) wave fronts. The asymptotic limit fre- 
quency of the tweek spectrum, however, is supposed 
to relate to horizontal wave fronts and should conse- 
quently be identical with our fundamental resonance 
frequency, at least in case of a well-defined refiection 
level in*^ the ionosphere. In other respects, the 
behavior of tweeks and that of our resonance 
phenomenon are not expected to be the sarne. 

The energy of the resonance oscillations is limited 
by losses. In our case the losses are due to partial 
transparency of the lower ionosphere (leakage), 
absorption on the path in the ionosphere below the 
reflection level, and imperfect ground reflection. 
Only the partial transparency is included in the 
theoretical treatment to follow. At the rather low 
frequencies the partial transparency is thought of as 
being the main cause of loss under most conditions. 
It was noted before that the ground is an almost 
perfect reflector. An estimate of the absorption 
in ionospheric propagation is made in appendix 1. 
The absorption coefficient found for normal daytime 
conditions may not prevent the resonance oscilla- 
tions, but we have to be aware that absorption tends 
to depress the quahty factor of the resonances and 
to smooth the resonance curves. This may prefer- 
ably happen in daytime and at the higher resonance 
frequencies. 

The refractive index of the ionosphere for our 
frequencies is in most heights fairly weU represented 
by the ''quasi-longitudinal" approximation (eq (14) 
of appendix 1, slightly simphfied) 



n'=lT 



n 



ffn I COS a I 



(1) 



The two signs refer to the ordinary ( — ) and the 
extraordinary (+) wave.^ /o is the plasma frequency, 
-fff the gyrofrequency of the electrons in the ter- 
restrial magnetic field, and a the angle between the 
wave normal and the direction of the terrestrial 
magnetic field. Equation (1), sometimes further 
simplified by neglecting unity, has frequently been 
used to describe whistler propagation [Storey, 1953], 
(more references given by Poeverlein [1959a]). The 
equation is based on the assumptions : /o>/, f^Jn, 
a not near it 12, and negligible colhsion frequency. 
In the lowest ionosphere there is, of course, a height 
interval in which /o is comparable to /. At our 
frequencies this height interval is found in a rather 
low altitude and is obviously narrow in terms of 
wavelengths (or decay distances). It has therefore 
no remarkable influence on the propagation and eq 
(1) may be used throughout. At higher frequencies, 
at which /o is comparable to / in a significant height 
interval, the refractive index formula has to be more 
rigorous and coupling or transition between ordinary 

1 This choice of signs complies with Storey's [1953] classical paper on whistlers 
and with the conventional definition of ordinary and extraordinary waves on 
basis of the Appleton-Hartree formula. The opposite notation, however, is 
occasionally found. 



and extraordinary waves has to be taken into con- 
sideration. This is known to be the case at fre- 
quencies of the order of 10 or 20 kc/s. 

Contributions of positive ions to the refractive 
index become noticeable near the gyrofrequency of 
the ions (approximately 50 c/s for 0+). All deriva- 
tions from eq (1) are therefore useful only for 
frequencies well above the ion gyrofrequency. 

1.3. Objective 

This theoretical study is concerned with the modi- 
fication of noise signals originating in the higher 
atmosphere by resonance of the space between earth 
and ionosphere. Frequency response curves of this 
space will be derived, wliich show for the frequencies 
under consideration the amplification or reduction 
of the signal intensity on the ground compared with 
the intensity of the incident wave. 

First some simplified models of the ionosphere are 
assumed, exhibiting one or two reflecting boundaries 
of homogeneous media. Later, samples of nighttime 
and daytime stratifications of the lower ionosphere 
are considered. With a simple refiecting boundary, 
resonance is obviously obtained when multiples of 
half-wavelengtlis fit into the resonance space. 

The study applies to wide horizontal wave fronts 
only. A noise emission coherent over a wide area is 
therefore required. Such an emission may arise from 
large-scale geophysical processes (e.g., incident cor- 
puscular radiation) rather than from local events. 

White noise or short pulses of a sufficient coherence 
over a wide area would lead to an immediate observa- 
tion of the derived frequency response curves. Mono- 
chromatic noise undergoes an enhancement or 
reduction of intensity depending on the position of 
the frequency in the frequency spectrum (response 
curve). Generally the noise spectrum observable on 
the ground is the result of the noise emission, eventual 
propagation effects on the path through the iono- 
sphere, and the response of the resonance space. 

Processes of noise generation in the higher at- 
mosphere have been treated to some extent in the 
literature [Gallet and Helliwell, 1959] (some more 
references are given by Gustafsson, Egeland, and 
Aarons [I960]). A survey of any amount of observa- 
tional material as for indications of our resonance 
effect has not yet been made. 

This study does not show much of a relationship 
with familiar theories of ELF or VLF propagation. 
This is mainly due to the consideration of waves 
that come down from high altitudes and travel ver- 
tically. Other investigations are concerned with 
propagation from a source to a receiver, both between 
earth and ionosphere. Some remarks on absorption 
and on a wave theory of propagation that might be 
of more general use appear in the appendixes. 

2. Simplified Models 

2.1. Earth and Sharply Bounded Ionosphere 

The earth will be considered as a perfect reflector 
in all cases that will be discussed. The simplest 
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model oi" the ionosphere is a lioiiiogeneous medium 
with a sliarp boundary at the bottom. Tlie refractive 
index of the ionosphere according to eq (1) has a high 
positive value for the extraordiiuiry wave and in 
general an imaginary value for the ordinary wave. 
The extraordinary wave will be dealt with now. 

The simple model of the resonance space and th(^ 
ionosphere is shown in figure 1. The electric liehl in 
the resonance space is drawn in for various resonance 
cases. The field may originate from some noise 
emission in the ionosphere. A descending wave 
generated in the ionosphere will reach the bottom 
of the ionosphere and lead to a standing wave in the 
resonance space below and to a reflected wave trav- 
eling back within the ionosphere. The standing wave 
in the resonance space may be interpreted as the 
superposition of a descending and an ascending wave 
traveling in a vertical direction. 

Ascending and descending waves will be denoted 
by F and G respectively. Strictly speaking, the 
symbols will denote the horizontal electric field 
strength components of these waves. In the aniso- 
tropic ionosphere an additional vertical (longitudinal) 
field strength component may exist. 

Figure 2 is a representation of tlie horizontal 
electric field strengths /''and G in the compkw plane. 
The vectors indicate amplitudes and j)hases of 
ascending waves {F) and descending waves {G). 
The polarization of tlie waves (disregarding vertical 
E components in the ionosphere) is assumed constant 
over the entire propagation path. In case of a 
reflecting artificial dielectric, for example, the polai- 
ization may be linear in a fixed dircM'tion. \\\ case 
of the extraordinary wave which is determined by 
the ionospheric propagation characteristics, the po- 
larization is circular in the quasi-longitudinal 
approximation that is being used now [Katcliffe, 
1959]. The sense of the rotation is clockwise for an 
observer looking in the direction of the terrestrial 
magnetic held. 

The construction of figure 2 starts at the earth's 
surface. The field strengths Fi and Gi there are 
opposite to each other iu consequence of the assumed 
perfect reflection. In traveling upward the two 
vectors rotate in opposite senses. The rotation in 
the air space amounts to kh and is, with a time factor 
exp {ioit), clockwise for the ascending wave {F), In 
traversing the bottom of the ionosphere there is a 
sudden change from F2 and G2 to F^ and G^. 
The total horizontal field strength components 

Ehor='F-\-G ^ 

H^or-=H=J'^n(F-G) f (^) 

\ Mo J 

must be continuous at the boundary, /t^or there is 
imaginary, i^hor real (see fig. 2). In passing the 
boundary, therefore, the end points of the vectors 
in figure 2 are shifted syjunietrically along a hori- 
zontal line with a reduction of the real compo- 
nents by a factor equal to the refractive index 
ratio (n of the ionosphere in our case). It may be 
recalled that n of the extraordinary wave is assumed 
to be real and greater than unity. 




E LAYER 



GROUND 



FiGUKE 1. Earth ami slxirply hounded ionosphere. 

The broken lines represent the distribution of the electric field strength in 
various resonance cases. The numbers denote the levels for which field strength 
vectors are shown in figure 2. 

It is superfluous to consider the longitudinal elec- 
tric field strength components of the waves at the 
boundary. Continuity of a I'esulting normal D com- 
ponent at any boun<hny is guarantcuMl hy requiring 
continuity of //,,()r- 

From oui- const ruction it is easily seen that the 
reduction of the field vectors F and G in going up 
tlirough the boundaiy is maximum if the vectors at 
th(^ bounih'U'X' ai'e r(\al or 



M^ 



:]7i 



(3) 



This is the resonance condition. In all cases given 
by e(i (o), a relatively strong field in the space 
between earth and ionosphere {Fi and (r,) ai'ises 
from a weak wave incident from above ((73). 

EcH'ause (73 is tfie field of the pi'imary wave, one 
might hk(^ to start the computation or construction 
of the field strengths with G^. This would howeve>r 
necessitate the consideration of multiph* refiections 
between earth and ionosphere, which did not become 
apparent in the method chosen. 

The procediu-e of figure 2 is also usabh' for more* 
complicated models of the ionosphere. A multi|)le 
layer with several boundaries rcniuires an it(M'ative 
application of the procedure (section 2.2 deals with 
a doubk^ layer). A step approximation allows us 
to transfer the method to a c;)ntinu()usly varying 
medium (section 3 and appendix 2). 

Figure 2, of course, is nothing but a graphical 
representation of the analytical expressions describ- 
ing the waves in the present case. From the figure 
or from the analytical fornnilation we derive the 
ratio between the field strcMigtlis arriving at the 
ground {Gi) and incident on the boundary of the 
ionosphere from above (6^3) 



IMAGINARY 




Figure 2. Electric field strengths in the complex plane (F 
electric field strength of the ascending wave, G electric field 
strength of the descending wave) . 
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■yjv} sin^ M+cos^ kh 



(4) 



The ratio of the Poynting vectors may be con- 
sidered as more significant than the ratio of the 
fiekl strengths. The horizontal field strength com- 
ponent G of the descending wave in a medium of 
refractive index n corresponds to a vertical com- 
ponent of the Poynting vector proportional to 
nG'^ (cf . eqs (2)) . An eventual horizontal component of 
the Poynting vector is not of interest in this study. 
The '^ratio of the Poynting vectors'' in the folio wing- 
always refers to the vertical components of the 
Poynting vectors which are proportional to nG^. 
The ratio of the Poynting vectors corresponding to 
the field strength ratio of eq (4) thus is 



n 



n^ sin^M+cos^M 



(5) 



Equations (4) and (5) show that resonance is obtained 
under the condition of eq (3), with an enhancement 
of the Poynting vector 



=n. 



(6) 



Maximum reduction of the Poynting vector is ob- 
tained for 

1 3 5 
2 



kh = ^ TT, ^ TT, 



TT, 



(7) 



witli a reduction ratio 



1 |^|2 

n \Gs\ " 



(8) 



For comparison the results of the Fresnel reflection 
theory for waves traversing a boundary may be 
noted: If the waves come from the medium with 
high-refractive index, a fraction 2n/(n+l) of the 
electric field strength or a fraction 4n/(7i+l)^ of the 
Poynting vector appears in the second medium. 
Equation (6) shows that the resonance effect in our 
case seems to improve the transparency of the 
boundary (in terms of Poynting vectors) by a factor 
of approximately n^j^. The reduction in the 
minimum case amounts to approximatel}^ 1/4 (eq 8). 

The enliancement of the Poynting vector in the 
resonance case does of course not mean a violation 
of tiie law of conservation of energy. With perfectly 
reflecting earth, there is no net power flux in the space 
between earth and ionosphere. The Poynting vector 
considered in this space is the Poynting vector of 
one of the two progressive waves of which the 
standing wave is supposed to be composed. 

2.2. Earth and Double-Layer Ionosphere 

The single-layer model (fig. 1) may be a fairly 
good approximation of the lower ionosphere at night. 




km 
GROUND 



Figure 3. Earth and double-layer ionosphere. 



The daytime ionosphere in lower heights is better 
approached by two adjacent dielectrics, representing 
the D and the E layers.^ This model is shown in 
figure 3. The heights are chosen in accordance with 
normal daytime conditions. The refractive indices 
for the extraordinary wave at a frequency around 
500 c/s can be assumed as 






= 50/ 



(9) 



The vector diagram, figure 2, must now be com- 
plemented by a rotation of the vectors corresponding 
to the phase shift in traveling up through the D 
layer and one more horizontal shift of the vector end 
points corresponding to the transition from the D to 
the E layer. In passing through the boundary 
between the D and E layers, the horizontal com- 
ponents of the field strength vectors in the vector 
diagram are reduced by the refractive index ratio 

UeIUj). 

In the present case two kinds of resonance may be 
distinguished, the one resulting from maximum 
reduction of the electric field strengths F and G 
(fig. 2) at the bottom of the D layer, the other 
from maximum reduction at the boundary between 
D and E layers. The first kind of resonance requires 
that the field vectors F and G in reaching the D 
layer are real, or in other words, that an integral 
number of half-wavelengths fits into the space 
between ground and D layer. With a height of 70 
km for the bottom of the D layer (fig. 3) this 
yields the resonance frequencies 



/-2.14X10^Xm(c/s) 
m=l, 2, 3, . . . 



}• 



(10) 



The minimums between the resonances are found at 
m=^, 2^, 2^, ... . 

The second kind of resonance, corresponding to 
maximum reduction at the boundary between the 
D and the E layer requires real field vectors at this 
boundary. With the model of figure 3 and the 
refractive indices of eqs (9) resonance frequencies of 
509, 1386, ... c/s and minimums at 151, 945, . . . 
c/s were derived. 

2 A brief account of a resonance case in the two-layer model was given in an 
earlier paper [Poeverlein, 1959b]. 
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A siiiii)l(' jipproxiniatioii for tJie second kind of 
resonance is obtained by assuming 7i/):^l. With 
large Ur, tlie liorizontal components of F and G (fig. 
2) in passing ttie bottom of the D layer are reduced 
to almost zero. This means that F and G imme- 
diately above the bottom of the D layer become 
imaginary. A phase shift of ^tt, ^tt, ^tt, . . . 
within tiie D layer makes then the vectors F and G 
real at the top of the D layer. The corresponding 
resonance frequencies under our conditions are 469, 
1406, 2340, . . . c/s. 

3. Realistic Cases 

3.1. E Layer, Extraordinary Wave 

A nighttime E layer profile given by Nertney 
([Poeverlein, 1959a] fig. 18) may now be used for the 
study of resonance of the extraordinary wave. The 
profile and an approximative step function are 
shown in figure 4. 

The refractive index in the ionosphere (eq 1) 
depends essentially on the plasma frequency /o. 
Values of /o following from 



fl=S0.7XlO^XN 



(11) 



are noted in figure 4 (/o in c/s, the electron concen- 
tration A^ per cm^). In the refractive index formula, 



/^|cosa| = 1.316X10^ c/s 



(12) 



may be assnmed (corresponding to some temperate 

latitude) and hence 



/oV(/H|cosa|)=61.3iV. 



(13) 



The step function of figure 4 similarly to a combi- 
nation of several dielectrics (sections 2.1 and 2.2) 
allows the iterative application of the graphical 
method (fig. 2). The step-by-step construction of 
t he field vectors was begun at the perfectly reflecting 
ground and conducted up to 104 km height. It 
appeared that not much contribution to the reflec- 
tion of the waves comes from heights beyond 104 
km, though some inaccuracy may arise from stopping 
at this somewhat low level and from the large width 
of the individual steps. 

fn,c/s 
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Figure 4. E layer profile. 
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FicuRE 5. Ualio of the Poynting vectors on the ground and in 
the ionosphere for the profile of figure 4- 

The graphical construction was used to derive the 
intensity variation of the descenduig wave. The 
ratio of the Poynting vectors on the ground and in 
104 km height, (l//^io4)X | 6^0/6^104! ^, is plotted vim'sus 
frequency in figure 5. Tlie positions of the maxi- 
mums and minimums are not much difFerent from 
those for a sharply bounded ionosphere (section 2.1), 
in 90 or 100 km height. The gradient of the electron 
concentration apparently is too steep to have much 
influence. A shallower gradient doubtless has some 
effect and must lead to diminished resonance fre- 
quencies as does a I) layer I)el()w the E reflection 
level (next section). The greatest enhancement of 
the Poynting vector (fig. 5) is seen to amount to 5 
(at 1.63 kc/s), the strongest [reduction to 0.11 (at 
0.60 kc/s). 
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Figure 6. Profile of D layer and E layer. 
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3.2. D Layer and E Layer, Extraordinary Wave 

In daytime, with presence of the D hiyer, the 
profile of tlie ionosphere may according to Pfister 
be that of fi^'ure 6 ([Poeverlein, 1959a], fig*. 17, 
curve (b)). The step approximation of this profile 
leads to the Poyuting vector ratio, (l/^ioo) X | 6^0/6^100 1 ^j 
represented in figure 7. The maximum enhance- 
ment is 1.02, tlie strongest reduction 0.06. 

The sequence of the resonance peaks seems to be 
irregular. In fact, the stratification of the iono- 
sphere resembles the simplified model of section 2.2 
for which two series of resonances were found. The 
first resonance peak in the present case (at 440 c/s) 
corresponds to the first peak of the simplified model 
(at 509 or 469 c/s) which was explained by the be- 
havior of the field at the bottom of the E layer and 
the phase shift in passing through the D layer. In 
the further course of the resonance curve (fig. 7) it 
seems impossible to separate the two series of reso- 
nance peaks of the simplified model. 

There is no pretension to accuracy of figure 7. 
At the higher frequencies, for which the subdivision 
of the step function (fig. 6) is not fine enough, some 
estimates were used in the values plotted. 

3.3. Ordinary Wave 

For the ordinary wave, which is obtained with the 
minus sign in eq (1), the refractive index is imaginary 
in most parts of the ionosphere. The graphical 
method of figure 2 is still usable for imaginar}" ri] 
however, the electric field strengths F and G become 
purely imaginary and have unequal amounts. The 
two waves {F and G) are attenuated, each in pro- 
gression in its direction of propagation. 
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Figure 7. Ratio of the Poynting vectors on the ground and in 
the ionosphere for the profile of figure 6. 



Tlie resonance case now is the case of disappear- 
ance of G in higher altitudes, strictly speaking in the 
altitude in which the signal is generated. In this 
case the enhancement of the Poynting vector is 
infinite, because no initial Poynting vector is present. 
A finite enhancement would result from considering 
electron collisions in the refractive index formula. 
The present study in whicli collision losses are dis- 
regarded allows therefore only the computation of 
the resonance frequencies. The first resonance fre- 
quency in the case without D layer (fig. 4) is found 
to be 1.52 kc/s, nearly the same as for a sharply 
bounded ionosphere. 

Wnereas resonance of extraordinary waves may 
bo excited by waves coming down from any altitude 
witiiin the ionosphere, the ordinary-wave resonance 
requires emission of radiation not too far beyond the 
reflection level, in order to have the ordinary waves 
only slightly attenuated in traveling down. The 
radiation source may be somewhat above the reflec- 
tion level of the ordinary waves or in the part of the 
ionosphere below the reflection level (lower tail of 
E, or D layer). The polarizations of ordinary and 
extraordinary waves show opposite senses of rotation. 

4. Comparison With Observations 

As noted in section 1.3, only an emission of white 
noise or short pulses in the ionosphere, fulfilling the 
coherence requirement, would lead to an observa- 
tion of the derived resonance functions (figs. 5 and 
7). In case of a given noise frequency the resonance 
function determines the enhancement or diminution 
of the signal on the ground compared with the noise 
intensity descending from higher altitudes. 

The field observable on the ground should be a hori- 
zontal FI. Oji a perfectly reflecting ground there 
should be no horizontal E. Imperfect ness of the 
ground reflection, however, makes a horizontal E 
possible. Circular polarization as expected from 
extraordinary (or ordinary) waves would correspond 
to two horizontal E (or H) components with 90-deg 
phase shift. 

Noise emissions exhibiting two or sometimes three 
peaks in the spectrum were observed in the auroral 
zone in connection with geomagnetic disturbances 
[Aarons, Gustafsson, and Egeland, 1960; Gustafsson, 
Egeland, and Aarons, I960]. The observed spectrum 
resembles at least roughly the resonance functions. 
The noise peaks in most cases were near 750 c/s and 
between 2 and 3.5 kc/s. The peak near 750 c/s is 
rather stable in frequency. In one event of noise 
emission, for example, it is reported to have vaiied 
as little as ±150 c/s. At other occasions there were, 
however, larger deviations. The observers suggest 
gyroresoiiance of protojis as cause of the 750-c/s peak. 

The second peak in the observed spectrum is more 
variable. It may be that this peak is caused by 
resonance enhancement of practically white noise 
and the first peak by gyroresonance radiation, pos- 
sibly in near coincidence with our first resonance 
peak under the prevailing conditions. A clue to the 
nature of the noises could be obtained from observa- 
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tioji ol I heir polarization. In proloji gyroradiation 
ordinary waves are emitted. Our standing-wave 
resonance, on the other hand, may preferrably be 
excited by extraordinary waves, as the remarks in 
section 3.3 indicate. Simultaneous noise observations 
on the gi'ound and in the ionosphere might also allow 
a decision about the character of the noises. 

For a discussion of VLF and ELF noises with 
I'cspect to characteristics and possible origin one may 
refer to Gustafsson, Egeland, and Aarons [I960]. 

5. Additional Remarks 

From a pulse-type signal emitted in the ionospliere, 
a series of pulses will be received on the ground. The 
first pulse reaches the ground immediately after 
leaving the ionosphere and traveling down through 
the air space, ^i'he others traveled back and forth 
between earth and ionosphere once or several times. 

As previously, oidy noise signals that liave exten- 
siv^o horizontal wave frojits are consich^'ed. The 
spectrum of t lie I'eceivcMJ seric^s of ])uls(»s is then neces- 
sarily identical with th(^ resonance s])ectruni obtained 
from white noise emissions in the ionosphei(\ be- 
cause a pulse (provided it is sliort) has the same spcu*- 
trum as white noise. In the simplified model of a 
sharply bounded homogeneous ionosphere with a 
refractive ijulex indepenckvnt of IVeciuency, there is 
no distortion of the pulse in tfu^ ionospheric reflec- 
tion, oidy a rcHluction of the field strength by a fac- 
tor {n — l)/{n-\-l). The geometric series of pulses 
thus obtained has. a Fouriei' sptM'truni correspojiding 
to eq (4) of the above resonance theory. This proves 
the equivalence of the two considerations hi this 
special case. 

The (juestion arises whether the whistlers observed 
in the audiofrequency band are also affected by 
nsonance of the space between (^ai'th and i()nos])here. 
Whistlers in coming down from liigh allituch^s may 
have extensive* wave fronts, but presumably not in a 
horizontal position. A possible variation of the 
resonance phenomena resulting from this inclination 
is not so much of interest now. The modification of 
the whistler by the resonance effect means nothing 
but a multiple repetition of the whistler signal, just 
the same as in the case of pulses. It will not be 
discussed now whether this repetition may or may 
not be recognizable in recordings. An important 
leature of tlie resonance spectrum is the first mini- 
numi, which appears in figure 5 at 600 c/s and will 
l)e slightly shifted towards higher frequencies in case 
of inclined wave fronts. This minimum may con- 
tribute to a low-frequency cutofiF of the observable 
whistler spectrum. 

The resonance effect this paper deals with refers 
to standing waves between two parallel reflectors. 
The curvature of the actual reflectors (earth and 
ionosphere) may necessitate some correction of the 
tlieory at the lowest frequencies.^ The lateral ex- 
tension of the field must be wide in terms of wave- 
lengths, but does not have to include the entire 
spherical shell. 

3 For the general theory of VLF propagation between curved earth and iono- 
sphere, refer to Wait [1960]. 



A resonance effect caused by the spherical shape 
of the earth is the resonance of the zero-order mode 
propagated between earth and ionosphere which results 
from winding-up of the waveguide to a spherical 
shell [Schumann, 1957; Konig, 1959; Wait, I960]. 
The fundamental frequency of this resonance effect 
corresponds to a phase variation of 27r along a path 
around the earth and amounts according to Schu- 
mann [1957] and Konig [1959] to approximately 
9 c/s. 

The fundamental frequency of our resonance effect 
in case of no leakage into the ionosphere is identical 
with the cutoff frequency of the first-order mode be- 
tween earth arul ionosphere, because this cutoff re- 
lates to verticahy traveling waves. The liigher-order 
resonance frequencies similarly represent the cutoff 
frequencies of the liigher-orchM" modes in this case. 

Appendix 1 . Effect of Collisions 

The quasi-longitudinal a|)i)r()ximation of the re- 
fractive iiuh^x (l(M'iv(Ml from the Api)leton-Hartree 
fornuda with iiu'lusion of the collision fi'equency v 
of the electrons is ^ [Katcliffe, 1959] 



m 



:1- 



co§ 



co(aj— ^Viajfl^lcos a.\) 



(14) 



This approximation is usahle at sufficiently high 
frequencies and at very low frequencies, provided 
the angle a between wave noi'mal and terrestrial 
magnetic field is not near 90 deg. In both cjises 
the approximative fornuila can be further simj)li(ied. 
The high-frequency approximation of n with I he 
assumption of co^coo, co;:$>co//, and v<^o) is 



2 co^ 



1 CO(V 



(15) 



At sufficiently low frequencies we may assume 
ojo »a),a;5>a;co^, j/<co^ and obtain then for the ex- 
traordinary wave 

,U f \V,_iA^^\ (16) 
Va'coiy |cos a|/ \ 2cOiy|cosa|/ 

The absorption in ionosphei'ic propagation is 
given by exp ( ^ Imag (n) • ds\ The absorption 
coefficient thus is at high frequencies 



"^ IT / \ I 

- I imag w| = 



1 COq ^ 

2c?'c 



and at sufficiently low frequencies 



(17) 



(18) 



The absorption coefficient for high frerjupiicios, eq 
(17), shows the famiUar frequency (U^pcndence of 

* Angular frequencies are written in place of the frequency quantities of eq (1). 



471 



the ^'nondeviative" absorption (with neglect of w^ 
compared with co). The absorptioD coefficient for 
very low frequencies, however, decreases with de- 
creasing frequency. This has the consequence that 
at our rather low frequencies the absorption in the 
lower ionosphere is only small. 

A numerical e\7aluation based on eq (14) and 
daytime data of the D layer (fig. 6 and v as noted) 
gave the following amounts of the absorption 
coefficient of the extraordinary wave at 1 kc/s (for 
/H|cosa|-1.2.10^ c/s): 

0.021 per km in 72 km with ^^=10^ per sec, 
0.019 per km in 80 km with j/=2-10^ per sec. 

These values should lead to a considerable absorp- 
tion on paths of 50 km length, but the absorption 
coefficient decreases rapidly above 80 km and the 
decisive vertical path lengths are shorter. 

A closer investigation of the absorption and its 
inffiience on the resonance curve are still desirable, 
but the above estimate suggests that tlie absorption 
is small and may not deteriorate badly the resonance 
curves, particularly not at the lower frequencies. 
At night the absorption coefficients are much smaller 
than the ones derived for da3'time and in daytime 
the lowest resonance frequency is below 1 kc/s, as 
figure 7 shows. 

It may be pointed out that the absorption co- 
efficients for high and very low frequencies are 
different not only in their variation with frequency 
but also in their dependence on cop. The absorption 
coefficient, of course, is not in an immediate connec- 
tion with the reflection coefficient of the lower 
ionosphere at very low frequencies. 

Appendix 2. Wave Theory Based on the 
Vector Method 

The construction of the complex field strength 
vectors (fig. 2) can be used for a rigorous wave-theo- 
retical treatment of propagation in a medium that 
varies continuously with one coordinate. It was 
seen that the method of constructing field vectors 
is applicable to a step approximation of a continu- 
ously varying medium. The start of the gi'aphical 
procedure was at a perfect reflector, in our case the 
ground. The solution consequently is a standing- 
wave solution. 

If propagation of plane waves in a continuously 
varying medium is to be studied, we replace the 
medium by an approximative step function and add 
a perfect reflector somewhere outside, thus creating 
conditions for the vector method and a standing-wave 
solution. A traveling-wave solution or any more gen- 
eral solution follows from superposition of two in- 
dependent standing- wave solutions. These two 
standing-wave solutions may either be obtained by 
assuming the perfect reflector at two different loca- 
tions or by assuming two alternative reflectors, one 
that does and one that does not reverse the electric 
field strength in reflection of the incident wave. The 
only approximation introduced is the step function. 



By making the steps differentials the method becomes 
rigorous. Nevertheless, the method has certain 
limitations. Before showing them, a brief analytical 
formulation will be given. 

The medium will be assumed to vary with the 
coordinate z. For waves propagated in the direction 
of the gradient, the field strengths may then be 
written 

£'^^(7(s)sin [\n{z)kd^+(t>{2)] 

Hy=iJ-n(z)C(z) cos [fn(z)hd^+ct>(z)] 

(19) 

Tlie time variation factor exp (io:t) has been omitted. 
The amplitude function C(z) and the phase function 
0(2) are sought. E^ and Hy represent transverse field 
strengths only. An additional longitudinal electric 
field component E^ may however exist, due to the 
anisotropy of the medium. The polarization of the 
wave in the x,y plane must be constant in order to 
have the continuity considerations of the next para- 
graph valid. The polarization may be linear (with 
Ex and Hy only) or circular (for example with 
Ey=iEx,Hx=-iHy). 

The continuity of £i and Hy at boundaries requires 
that at the steps of the step function 



A(Cs) = 
MnCc) = 



}■ 



(20) 



Herein s and c are abbreviations for the sine and 
cosine appearing in eqs (19) : 



s=sin [ nkQdz-\-(l)] 
c=cos [ nhdz+(l)] 



(21) 



The integral term in tlie argument of these functions 
varies continuously. Only </> is supposed to jump at 
the steps of the step function. The differences of 
-s' and c needed in eqs (20) are therefore 



} 



As=c/\(l) 
Ac=—sA(l) 

Equations (20) together with eqs (22) yield 

AC An 



(22) 



c 



n 



A6= — sc 
n 



(23) 



Making the steps differentials leads to the rigorous 
wave-theoretical differential equations for C and (p 



dC^_dJi 

C n 

d6= — sc 
n 



(24) 
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It has to be noticed that the quantities s and c 
in eqs (24) are functions of </> (eqs (21)). The last 
of (Hjs (24) allows the coin])utation of the phase 
rmiction without considering- tlu^ amplitude func- 
tion, whereas the first equation involves both the 
ani])litu(le function and the phase function. Equa- 
tions (24) are also readily obtained from Maxwell's 
etjuations. 

A set of equations closelv related to eqs (24) has 
been derived by van Cittert [1939] [Kofink, 1947]. 
Van Cittert considered oblique incidence on a one- 
dimensionally stratified nunhuin. His equations 
therefore are more general in this respect. The 
medium he assumes as isotropic. The equations are 
necessarily different for the two cases of (linear) 
polarization. 

Slow variation of the refractive index suggests av- 
eraging the factors c^ and sc over many periods of 
c and .s. Equations (24) then become approximately 



d/l__ldn'\ 
C 2 






(25) 



wlience wc obtain 



(7 ex 4= 

■yjn 

<^= const 



(26) 



This is the behavior of amplitude and phase in tlic 
geometric-optical approximation. 

The above study was limited to nonnal incidence 
on a one-dimensionally stratified medium . A general- 
ization for oblique incidence (on a one-dimensionally 
stratified medium) may not be difficult, as van 
Cittert's formulation shows. TIk^ more severe limita- 
tion arises fi'om the demand of constant polarization. 
In case of oblique incidence the polarization in the 
planes of the stratification has to stay constant. 
This is easily achieved in isotropic media. In 
anisotropic media constant polarization is obtained 
in normal incidence if the polarization can be as- 
sumed as circular, independent of the location in the 
layer. In more general cases of propagation in 
anisotropic media (oblique incidence or elliptic 
polarization) the polarization in the planes of the 
stratification varies. The continuity at the steps of 
our step function can then be preserved only by 
simultaneous consideration of ordinary and extra- 
ordinary waves witli variable amplituck^ factors. 
This indicates the well-known phenomenon of cou- 
pling between the tw^o types of waves. 



The step method is not apphcable near a zero oj the 
refractive index. In approaching the zero, the re- 
quired step length decreases towards zero. However, 
a rigorous expr(»ssion for the refractive index that 
includes the collision frequency has no zero at any 
location. 

Complex refractive index leaves the analytical 
formulations valid, but makes the graphical pro- 
cedure of figure 2 more intricate. 



The author thanks James R. Wait for a stimulating 
discussion and Thomas J. Birmingham for help in 
preparing the manuscript. 
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